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LOCAL SYNTHESIS AND TOOTH CONTACT 
ANALYSIS OF FACE-MILLED, UNIFORM TOOTH 
HEIGHT SPIRAL BEVEL GEARS 


by 

F. L. Litvin 1 and A. G. Wang 2 
Department of Mechanical Engineering 
University of Illinois at Chicago 
Chicago, IL 


ABSTRACT 

Face-milled spiral bevel gears with uniform tooth height are considered. An approach 
is proposed for the design of low-noise and localized bearing contact of such gears. The 
approach is based on the mismatch of contacting surfaces and permits two types of bearing 
contact either directed longitudinally or across the surface to be obtained. Conditions to 
avoid undercutting were determined. A Tooth Contact Analysis (TCA) was developed. This 
analysis was used to determine the influence of misalignment on meshing and contact of the 
spiral bevel gears. A numerical example that illustrates the theory developed is provided. 
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NOMENCLATURE 


o. g Blade angle of gear head cutter (fig. 4) (Table 2) 

a v Profile angle of pinion head cutter (figs. 7, 18)(Table 3) 

7 i >72 Angles of pinion and gear pitch cone, respectively (figs. 6, 11, 12) (Table 

1) 

7 Shaft angle (Table 1) 

Tj\ Tangent to the path of contact on the pinion surface (Table 3) 

0 P Surface parameter of the pinion head cutter 

0 g Surface parameter of the gear head cutter 

A p Surface parameter of the pinion head cutter (figs. 7, 18) 

0 i 2 Angle formed between principal direction ej and e s (fig. 10) 

<j>i (z = 1,2) Angle of rotation of the pinion (z = 1) or gear (i = 2) in the process of 
meshing (figs. 11, 13, 14, 15) 

V'c, (* = 1, 2) Angle of rotation of the cradle in the process for generation of the pinion 
(t = 1) or gear (z = 2) (fig. 5) 

r/>,- ( i = 1,2) Angle of rotation of the pinion (i = 1) or gear (i = 2) in the process 
for generation (figs. 6, 8) 

u>i (z = 1,2) Angular velocity of the pinion (z = 1) or gear (z = 2) (in meshing and 
generation) 

w Ci (z = 1,2) Angular velocity of the cradle for the generation of the pinion (z = 1) 
or gear (z = 2) 

£, (z = 1,2) Pinion (z = 1) or gear (z = 2) tooth surface (fig. 17) 

S (i (z = 1,2) Pinion (z = 1) or gear (z = 2) generating surface (figs. 1, 2, 3, 16) 

AA P , A A g Pinion and gear axial displacements, respectively (figs. 11, 13)(Table 

4 ) 

AE, Ay Errors the offset and shaft angle, respectively (figs. 12, 13)(Table 4) 


v 



A<fo(<£ 1) 




h d 

kf , kh. , i ^q 

L,. 

m 21 

M 

Ma 

n fc> , N fc . 


iV t - (i = 1,2) 
9 , (*' = 1 , 2 ) 




Function of transmission errors (figs. 14, 15) 

Unit vectors of principal directions of pinion and gear tooth surface, 
respectively (fig. 10) 

Dedendum height of the pinion 

Principal curvatures of the pinion and gear tooth surfaces, respectively 
Matrix of orientation transformation from system Si to system Sj (3X3) 
Derivative of (Table 3) 

Mean contact point (figs. 1, 2, 3, 7, 9)(Table 3) 

Matrix of coordinate transformation from system Si to system Sj (4X4) 
Unit normal and normal to the generating surface £,■ represented in 
coordinate system Sk 

Number of teeth of pinion (i = 1) and gear ( i = 2) (Table 1) 
Installment angle for the head cutter of the pinion (i = 1) and gear 
(i = 2) (fig. 5)(Tables 2, 3) 

Radius of the generating surface of revolution for the pinion (figs. 3, 7, 


16, 18)(Table 3) 

Rj,, R g Radius of the head-cutter at mean point for the pinion and gear (figs. 

1, 2, 3, 4, 7, 16, 18)(Tables 2, 3) 


r t - Position vector in system 5, (i = 1, 2, h, < 1 , < 2 ) 

S Ti (i =1,2) Radial setting of the head cutter of the pinion (i = 1) and gear (i = 2) 
(fig. 5) (Tables 2, 3) 

Si Coordinate system 

vM (1 = 1,2) Velocity of contact point in its motion over surface £ t| 

vW) Relative velocity at contact point (i,j = 1,2, Ci,C 2 ,fi,t 2 ) 

up), Components of the velocity of the contact point in its motion over £1 


Vl 



1 Introduction 


Two models for spiral bevel gears with uniform tooth height were proposed by Litvin et al. 
[1]. The generation of tooth surfaces of such gears is based on application: (i) of two cones 
that are in tangency along their common generatrix (model 1), and (ii) a cone and a surface 
of revolution that are in tangency along a common circle (model 2). The pinion and the 
gear are face-milled by head-cutters whose blades by rotation form the generating surfaces. 

The generating surfaces provide conjugate pinion-gear tooth surfaces with a localized 
bearing contact that is formed by a set of instantaneous contact ellipses. The path of 
contact is directed across the surfaces in model 1 (fig. 1), and in the longitudinal direction 
in model 2 (fig. 2). The transmission errors are zero but only for aligned gear drives. 

It is well known that misalignment of a gear drive causes a shift of the bearing contact 
and transmission errors. The transmission errors are one of the main sources of vibration. 
Therefore, the direct application of the models discussed above for generating surfaces is 
undesirable. 

It was discovered that misalignment of a gear drive causes an almost linear but discon- 
tinuous transmission function. However, such functions can be absorbed by a predesigned 
parabolic function of transmission errors. The interaction of the parabolic function and a 
linear function results a parabolic function with the same parabola coefficient [2]. Based on 
this consideration, it becomes necessary to modify the process discussed above for generation 
to obtain a predesigned parabolic function of transmission errors. It was proposed in the 
work [3] to obtain the desired parabolic function of transmission errors by executing proper 
nonlinear relations between the motions of the cradle and the gear (or the pinion) being 
generated. This approach requires the application of the CNC machines. 

The purpose of this report is to propose modifications of generating surfaces that will 
obtain: (i) a localized bearing contact that may be directed in the longitudinal direction 
or across the surface, and (ii) a predesigned parabolic function. These goals, that will be 
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proven later, are obtained by the proper mismatch of the ideal generating surfaces shown in 
figs. 1 and 2. The mismatch of surfaces is achieved by application of modified generating 
surfaces shown in fig. 3. The modified generating surfaces are in point contact instead of 
tangency along a line that the ideal generating surfaces have. The desired parabolic function 
of transmission errors, the orientation of the path of contact, and the magnitude of the major 
axis of the contact ellipses are obtained by the proper determination of the curvature and 
the mean radius of the surface of revolution of the generating tool. 

Design of drives with a small number of pinion teeth may be accompanied with pinion 
undercutting. Using the approach proposed in [7, 8, 9], it becomes possible to avoid under- 
cutting of spiral bevel pinions. The meshing and contact of the tooth surfaces was simulated 
by the TCA (Tooth Contact Analysis) computer program developed by the authors. 

The contents of the report cover the following topics: 

(1) Method for generation of conjugate pinion-gear tooth surfaces. 

(2) Derivation of gear and pinion tooth surfaces. 

(3) Local synthesis as the tool for the directed mismatch of contacting surfaces. 

(4) Simulation of meshing and contact of misaligned drives. 

(5) Avoidance of pinion undercutting. 

Numerical examples for the illustration of the proposed approach are considered. 
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2 Method for Generation of Conjugated Pinion- Gear 
Tooth Surfaces 


Gear Generation: 

The head-cutter for gear generation is provided with inner and outer straight-line blades 
(fig. 4), that form two cones while the blades are rotated about the Zt 2 - axis of the head 
cutter. These cones will generate the convex and concave sides of the gear profile, respectively 


[12]- 


We apply coordinate systems S C2 , S 2 , S m that are rigidly connected to the cradle of 
the generating machine, the gear and the cutting machine, respectively (figs. 5 and 6 ). The 
cradle with coordinate system S C2 performs rotation about the Z m -<ixis : and ip C2 is the current 
angle of rotation of the cradle (We take i = 2 in the designations of fig. 5 ). Coordinate 
system St 2 is rigidly connected to the gear head-cutter that is mounted on the cradle. The 
installment of the head-cutter is determined with angle q 2 and S r2 = \0 C2 O t2 (fig. 5 (b)). 
The gear in the process for generation performs rotation about the Z b - axis of the auxiliary 
fixed coordinate system S b that is rigidly connected to the S m coordinate system (fig. 6 ). 
The installment of S b with respect to S m is determined with angle 72 , where 72 is the angle 
of the gear pitch cone. The current angle of gear rotation is ip 2 (fig. 6 ). Angles ip C2 and t /> 2 
are related as 

~r - = sin 72 ( 1 ) 

y>2 w 2 K ’ 

The observation of this equation guaranties that the AT m -axis is the instantaneous axis 
of rotation of the gear in its relative motion with respect to the cradle. 

Pinion Generation: 


The head-cutters for pinion generation are provided with separate blades that will gen- 
erate the convex and concave sides of the pinion profile, respectively (fig. 7 ). The pinion 
generating tool is installed on the cradle similarly to the installment of the gear generating 
cone (We take i = 1 in the designations of fig. 5). An auxiliary fixed coordinate system 
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S a is rigidly connected to the S m coordinate system (fig. 8). An imaginary process for the 
pinion generation for the purpose of simplification of the TCA program is considered. The 
installment of coordinate system S a with respect to S m is determined in the real process of 
cutting by the angle 71 that is measured clockwise, opposite to the direction shown in fig. 
8. The pinion performs rotation about the Z a -axis and ipi is the current angle of rotation. 
The angles of rotation of the pinion and the cradle are related as 


0C1 


— - = sm 71 


( 2 ) 


Axis X m in accordance to equation (2) is the instantaneous axis of rotation of the pinion 
in its relative motion with respect to the cradle. 
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3 Derivation of Gear Tooth Surface 


Equations of Gear Generating Surface 

We consider that the gear head-cutter surface is represented in St 2 (fig. 4) by vector 
function r t2 (s g ,9 g ) 

( R g - s g sin a g ) cos 9 g ' 

= (R a sin ot g ) sin 9 g (3) 

Sg COS OLg 

where s g and 9 g are the surface coordinates; a g is the blade angle; R g is the radius of the 
head-cutter at the mean point. Equations (3) may also represent the convex side of the 
generating cone considering that a g is negative. 

Coordinate system S t2 is rigidly connected to coordinate system S C2 , and the unit normal 


to the gear generating surface is represented by the equations 


n ($ ) — _^£2_ n _ §£h x 
cz( 3) |N C2 | ’ C2 d9 g ds g 


Equations (3) and (4) yield 


cos a g cos 6 g 
n C7 (9 g ) = cosa g sin9 g 
sina p 


(4) 

(5) 


Equations of the Family of Generating Surfaces in S2 

A family of tool surfaces is generated in gear coordinate system S 2 while the cradle and 
the mounted tool and the gear perform the rotational motions that are shown in figs. 5 and 


6. The family of surfaces is represented in S 2 by the matrix equation 

9 g , ip^) = M2fc(02)M6 m M m c 2 (^’c2)Mc2t2r(2(^? 9g) , . 

= (b) 

The product of matrices M 2 t 2 is based on the coordinate transformations from S t2 to ^2 

(figs. 5 and 6), where and ip C2 are related by equation (1) and 


M 26 = 


cos r />2 
— sin xp 2 
0 
0 


sin \j >2 0 0 

cos V> 2 0 0 
0 1 0 
0 0 1 


(7) 
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h^bm — 


sin 72 0 —cos 72 0 

0 10 0 

cos 72 0 sin 72 0 

0 0 0 1 


Mmc2 — 


cos V> C2 
sin V> C2 
0 
0 


— sin ip C2 0 0 

cos t/'c 2 0 0 

0 1 0 

0 0 1 


^^C2t2 — 


10 0 S r2 cos q 2 
0 10 5r2Sin92 
0 0 1 0 
0 0 0 1 


Equation of Meshing 

We derive the equation of meshing between the generating surface and gear as 


n c2 • v<? 2) = f(s g , 8 g , V> 2 ) = 0 


representing the vectors in S c 


where v^ j2 ^ is the relative velocity that is represented in the coordinate system S CJ ■ Here, 

V C2 22) = ui c? 2) x r c2 = (^C2 2> “ u c?) x r CJ (12 


where 




= M C 2 < 2 r t2 


u> i sin 71 


— s s sin Q g cos 0 g + B\ 
—s g sin a g sin 6 g + B 2 

Sg COS Ctg 


0 1 r — Ni/N 2 COS ~f 2 cos iJj C7 
-L Cjm L mb L b2 0 = Ni/N 2 cos 72 sin t/> C2 

u > 2 0 


Using the designations 


Bi = R g cos 6 g + S r2 cos q 2 
B 2 = R g sin d g + S r2 sin q 2 


and considering that N x /N 2 = sin ji / sin 72 and |wi| = 1 , we obtain from equation( 12 ) that 

s g cos 72 sin t/> C 2 cos a g N x / N 2 

= s s cos cos cos a 9 ^\ /N 2 ( 16 ) 

_ Ni/N 2 cos 72 {s g sin a g sin ( 8 g + ip c 2 ) — B x sin - B 2 cos xp C2 ) _ 
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The equation of meshing (11) is represented as 


s g{@gi ^ ci ) 


sin ctg(Bi sin 0^ + j? 2 cos 0 C2 ) 
sin(0 s + 0^) 


(17) 


Equations of Gear Tooth Surface 

Equations (6) and (17) represent the gear tooth surface by three related parameters. 
Taking into account that these equations are linear with respect to s g , we may eliminate s g 
and represent the gear tooth surface by two independent parameters, 6 g and 0 2 , as 


r 2 = r 2 (0 s ,0 2 ) 


(18) 
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4 Derivation of Pinion Tooth Surface 


Equations of Pinion Generating Surface 

The derivations are similar to those that have been described in section 3. The generating 
surface of revolution is represented in S tl (fig. 7) as 


r *i (^j>* &p) ~ 


[Rp - i?i(cosa p - cos(a p + A p ))]cos0 p 
[Rp - ili (cos a p — cos(a p + A p ))] sin0 p 
-ill (sin a p — sin(a p + A p )) 


(19) 


where A p and 6 V are the generating surface coordinates; a p is the profile angle at point M; Rp 
is the radius of the head-cutter at mean point; ill is the radius of the surface of revolution. 
Equations (19) can also represent the concave side of the generating surface of revolution if 


we substitute a p as 180° — a p . 


Coordinate system St 1 is rigidly connected to coordinate system 5 Cl , and the unit normal 
to the pinion generating surface is represented by the equations 

N c 


n Cl (A p ,0 p ) = j^L 


_ dr tl dr h 

C1 d0 v d\ , 


(20) 


Equations (19) and (20) yield 


n Cl (A p , 9p) — 


cos Op cos(a p + A p ) 
sin Op cos(a p + A p ) 
sin(a p + A p ) 


( 21 ) 


Equations of the Family of Generating Surface 

The family of generating surfaces is represented in Si by the matrix equation 


Fl(A p , Op, t/^l) — ^^la(0l)h^ a mh^[ m ci (V^ci (A p , 0 p ) 

= Mitj (A p , 9p) 


( 22 ) 


The product of matrices Mi tl is based on the coordinate transformations from 5 ( , to Si 
(figs. 5 and 8), where tpi and ip Cl are related by equation (2). Here, 


Mia 


cos -01 — sin 0i 0 0 
sin ipi cos 0i 0 0 

0 0 10 

0 0 0 1 


(23) 
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— sin 71 0 — cos 71 0 

0 10 0 

COS71 0 — sin7i 0 

0 0 0 1 

cos V'cj — sin V’c, 0 0 

sint/> Cl cos^ci 0 0 

0 0 10 

0 0 0 1 


M C1 tj 

Equation of Meshing 

We derive the equation of meshing between the generating surface and pinion as 

n ci ' v cj l1 ^ = /(A p ,0p, t/>i) = 0 (27) 

where v^* 1 ) is the relative velocity in the coordinate system S Cl . The vectors are represented 
in S Cl . Here, 

vg 11 ’ = <■’ X r t , = «>-«<,") X i*> (28) 

where 

r Cl =M Cltj r ll 

Ri cos(q p + A p ) cos 0 P + ] 

= Ri cos(q p + A p ) sin 0 P + B 2 ' ' 

— i?i(sina p - sin(ap + A p )) 

0 0 u)\ cos 71 cos ip Cl 

(<*>€? } -w l 1) ) = 0 - L Cim L ma L al 0 = —<jJi cos 7 x sin ip Cl (30) 

u>\ sin 7 a J [ — J 0 

Using the designations 

B\ = Rp cos 0 P + S T 1 cos q\ — R\ cos a p cos 6 P 
B 2 = Rp sin 0 P + S rl sin qi — R x cos ct p sin 0 P 

and considering that |u>i| = 1 , we obtain from equation(28) that 

R\ (sin q p — sin(a p + A p )) cos 7 j sin ip Ci 
v^ lJ) = i?i(sina p — sin(a p + A p ))cos 7 ! cost p Cl (32) 

. cos 71 (Ri cos (ctp + A p ) sin(0 p + t p Cl ) + £1 sin xp Cl + B 2 cos ip Cl ) _ 




(24) 

(25) 
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The equation of meshing (27) is represented as 


tan( q p + A p ) 


—Ri sin q p sin(fl p -f V> Cl ) 
B\ sin tp Cl + B 2 cos V’d 


(33) 


Equations of Pinion Tooth Surface 

Equations (22) and (33) represent the pinion tooth surface by three related parame- 
ters. After elimination of parameter A p we may represent the pinion tooth surface by two 
independent parameters, 0 P and ii>\. 


ri =r!(0p,t/>i) 


(34) 
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5 Local Synthesis 

The ideas of local synthesis are based on the following considerations [2]: 

(1) The pinion and gear tooth surfaces are in tangency at the mean contact point M that 
is in the middle of the contacting surface. 

(2) The gear ratio is equal to the theoretical one. 

(3) We have to provide in the neighborhood of M the following transmission function 

(fig- 9) 

<W.) = (35) 

where \rn' 2 i is the parabola parameter of the predesigned parabolic function of transmission 
errors 

Afcfo) = (36) 

(4) In addition it is necessary to provide the desired direction of the contact path. 

All these goals can be achieved by the proper mismatch of the contacting surfaces of the 
pinion-gear tooth surfaces. The solution to this problem requires directions of the contacting 
surfaces. However, since the equations of the pinion and gear tooth surfaces are represented 
in a complex form, we will represent the principal curvatures and directions of the generated 
surfaces in terms of the principal curvatures and directions of the generating surfaces (the 
head-cutter surfaces) and the parameters of motion. The procedure of the local synthesis is 
as follows: 

Step 1: We consider as given the surface of the head-cutter that generates the gear 
tooth surface. The head-cutter surface is a cone and is in line contact with the surface of the 
gear. One of such contact lines passes through the mean point M of tangency of the pinion 
and the gear tooth surfaces. Considering the surface of the gear head-cutter as known, we 
determine at point M the principal curvatures and directions of the gear head-cutter. 
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The principal directions of the generating cone are 

e / t2) = tr/tr! = (- sin ^ cos ^ °l r 


(ta) r /j • • /) it 

e^ 2 ' = ~ I - sin a 9 cos -sinajSinfl, cosQjJ- 1 

The principal curvatures of the generating cone are 

k = — cosa<,/(i?<, — s ? sina 5 ) 

( 38 ) 

*l w = 0 

Step 2: Our next goal is to determine at M the principal curvature k„ and k q and the 
principal directions of the gear tooth surface £2- We apply for this purpose the equations that 
have been proposed in [2] and represent the direct relations between the principal curvatures 
and directions for two surfaces being in line contact. 

Surface £ tj and £2 are in line contact when cone £ t2 generate the gear tooth surface £2. 
The principal curvatures of the gear k s and k q can be obtained from the equations 

ta ° 2<T " = 


tan 2a g = 


k, + k, = + k\'_" + &L+J1 


where 


6,3 = -kf ) vf 2) + [nw<‘»»>ey ,) ] 

6*, = + [nu^’M’ 1 ] 

633 = -kf\vf 2) f - + [nwO^vO* 2 )] — n-[(ui^ 2 ^ x V(^) — (u>( 2 ) x vj* 2 ^)] 

( 40 ) 

The principal directions on the gear tooth surface are represented by unit vectors e, and 


e„, where 


s _ cos a g sin a g 
q — sin a g cos a g 
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Step 3: We now consider that the gear and pinion tooth surfaces, E 2 and Ex, are in 
tangency at M. As a reminder, the mismatched gear and pinion tooth surfaces are in point 
contact at every instant. 

Unit vectors e, and e, represent the known directions of the principal directions on 
surface E 2 . The principal curvatures k s and k q on the gear principal directions are known. 
Our goal is to determine angle a 12 that is formed by vectors e/ and e, (fig. 10) and the 
principal curvatures kf and kh of the pinion tooth surface at point M. Unit vectors e/ and 
represent the sought-for principal directions on the pinion tooth surface Ex- 
Step 4 : The three unknowns: kj, kh and <r 12 can be determined using the approach 
developed in [2]. We use for this purpose the following system of three linear equations. 
Three linear equations that related the velocity v(U of the contact point over surface Ex are 
derived in reference [2] as: 


a n v 
a i2V. 


P + 
« + 


a 12 v 


(i) = 

h = 


0X3 

023 


(42) 


+ a>2 3Vy> = <233 


-J- n-ioiJ 1 ) = 

s ' 

The augmented matrix formed by the coefficients a,x, a l2 and a t3 is a symmetric one [2], 
Here, an d v (i) are the components of the velocity of the contact point that moves in 
the process of meshing over the pinion tooth surfaces Ex- Coefficients a tl , a l2 and a l3 are 
represented in terms of k 3 , k q , k j , kh, cr 12 and the parameters of motion. 

Step 5: Equation system (42) represents a system of three linear equations in two 
unknowns: and . Surface Ex and £ 2 are in point contact, the path of contact has 

a definite direction, and the solution of equation system (42) with respect to and 
must be unique. Therefore, the rank of the augmented matrix formed by a tl , a, 2 and a, 3 is 
equal to two. This yields that 


0ll022 a 33 + Ul2 a 23°13 + Ul3< I 12 a 23 — a 22 a l3 ~ a ll a 23 ~ a 33 a x 2 


flu 

Ux2 

^13 

ai2 

a 22 

a 23 

^13 

& 23 

a 33 


= F{kf, kh, k s , k q , (T 12 , m 21 ) = 0 


(43) 
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Here, 


a u = k s — kf cos 2 &12 — khsin 2 cri2 
a 12 = 0 .h(kf - fc/ l )sin 2 cTi 2 
ai3 = -ksvW + [nw ( 12 ) e a ] 

a .22 = — fc/ sin 2 (Ti 2 — ^ cos 2 (Ti 2 (44) 

a 2 3 = -M? 12) + [nu;( 12 )e,] 

033 = M W S I2) ) 8 + ^?(^ 12) ) 2 ~ [ nu ?( 12 W 12 >] 

— n KwW x vg>) - (u/ 2 > x v t (2) )] + m' 21 ( n x k 2 ) • r m 

a 12 — cr g — a p is the angle between the principal directions of this two contacting surface, 
m' 21 is the derivative of <fo(<M at the contact point. Coefficient 033 contains the derivative 

m! ' 21 = («) 


where 


m 2 1 = 


d^2 

d<f>\ 


(46) 


From equation (43), we get 


033 = 


— Q12Q23Q13 ~ Q13Q12Q23 + Q22O13 + Qll Q 23 
O11O22 _ o 2 2 


(47) 


Substituting equation (47) into equation (44), we can obtain the derivative 

033 - kM 12) Y - kM 12) ) 2 + M 12 M 12 )] + n-[{ U W X v t (2) ) - (u/ 2 ) x v t ( r 1} )] 


m 21 = 


(n x k 2 ) • r„ 


(48) 


The parabola coefficient of the parabolic function (36) is m! 2x j 2. 

The other relation between the coefficients o tl , a t2 and a l3 may be determined considering 
that 

t,(i) 

tan»7i = (49) 

vi ’ 

where r/i is the assigned direction at M of the tangent to the path of contact on the pinion 
surface Ej. 

Using the relations discussed above between the coefficients of linear equation [2], we are 
able to determine the sought-for pinion principal curvatures kf,kh and orientation angle <Ti 2 . 
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Step 6: We consider now that for surface Ei the followings are known: (i) the principal 
directions determined by unit vectors e/ and (ii) the principal curvatures kj and kh, 
and (iii) angle <Ti 2 formed by unit vectors e/ and e 3 (fig. 10). Our goal is to determine 
the principal curvatures and directions of the pinion head-cutter generating surface that 
is designed as the surface of revolution (fig. 3 ). The pinion head-cutter surface and the 
pinion tooth surface are in line contact at every instant. Using the direct relations between 
the principal curvatures and directions for two surfaces being in line contact [2], we may 
determine the principal curvatures and principal directions of the pinion head-cutter. Then, 
the desired mismatch of the surfaces of the gear and the pinion will be provided by the 
generation of the gear and the pinion by the designed head-cutters. 

Using the approach discussed above, we obtain the following equations 
—2b\r,b'o-, 


tan 2<7 P 

V 2 U 2 


°23 ~~ °13 ' 

kh ~ k f 

— -2613623 

632 sin 2^ 

kh + kf 

= kw + kp 


*33 


( 50 ) 


^33 


where 


b'13 = — 

*23 = -k^vM + [nw<*i>e<«»>] 


b ' 33 = -k^)(v^y - *<**>( t;^ 1 )) 2 + [nu^M 111 )] - n-[(u/<*i> x vg>) - («(») x vg 1 *)] 

( 51 ) 

The principal directions on the pinion and the pinion head-cutter are related as follows 


e / 


COS (Tp 

sin a p 


e (h) 1 

5 

. e/i . 


— sin cr p 

COS (Tp 


e (‘i) 

J 


( 52 ) 
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where the principal directions of the generating surface of revolution are 
e > ,,) = = ! - sin "' c ° s ^ ° f 

e<*‘> = = [- sin(a p + A p ) cos & v - sin(a p + A p ) sin0 p cos(a p + A P )] T 

<7A p (/Ap 

(53) 

The principal curvatures of the generating surface of revolution are 

k M = - cos(q p + X P )/(R P + i?i(cos(a p + A p ) — cos a p )) 

(54) 

*? l) = -l/Rx 

Equations (53) and (54) permit the representation of the principal curvatures and direc- 
tions on the pinion head-cutter surface in terms of R\ and Rp. We remind that equations 
(48) and (49) contain parameters R x and Rp (figs. 2). Considering as given m' 2x and 771, we 
can determine from equation (48) and (49) R\ and Rp. 

Step 7: At this step we know the principal curvatures and directions on the pinion and 
the gear, and the principal curvatures and directions on the pinion and gear head-cutters. 
The obtained mismatch of pinion and gear tooth surfaces will provide in the neighborhood of 
the mean contact point the desired parabolic function of transmission errors and the direction 
of the contact path. The principal curvatures and directions obtained on the pinion and 
gear head-cutters will provide the required mismatch of the pinion and gear tooth surface. 
Our next goal is to determine the dimensions of the instantaneous contact ellipse and its 
orientation, considering as given the elastic approach of the contacting surfaces. The solution 
is based on the following procedure [2]. 

The major axis and minor axis of the contact ellipse can be determined as 



where S is the elastic approach obtained from experimental data; A and B are determined 
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by 

A = ~ k ¥ } - (9i - 2 9 i92 cos2<7!2 + gl)*) 

B = 4l fc £ 1) - + ( 9i - ^9\92 cos2<7 12 + gj)^} 

and 

k z ] = kf + k h k £ 2) = k a + kg 

9i - kf - k h g 2 = k 3 - k q 

The orientation of the contact ellipse in the tangent plane is determined by 


(56) 


(57) 


cos 2a 


<7l — 92 COS (7i2 


(58) 


(^i ~ ^9\92 cos 2 oi 2 + gl)'* 

Directions for the Computational Procedure of Local Synthesis 

Step 1: The parameters of the gear head-cutter and its installment are considered as 
known (see, for instance, Table 2). 

Step 2: The mean contact point is considered as known as well (It is determined by the 
application of the TCA program that provides the tangency of contacting surfaces at the 
mean contact point). 

Step 3: Using equations (39) and (41), we determine at the mean contact point the 
principal curvatures (k 3 , k q ) and the principal directions represented by unit vectors e, and 
e„ (fig. 10). 

Step 4: We use the values Ri and R v (fig. 3) as the first guess for the pinion head 
cutter. Then, we determine at the mean contact point the principal curvatures k^ a l \ k^ ll \ 
and principal directions represented by e^ 1 ^, e^d applying for this purpose equation (53) 
and (54). 


Step 5: We compute the principal curvatures kf, kh of the pinion tooth surface and 
principal directions represented by unit vectors e / and e/, applying for this purpose equations 
(50) and (52). 

Step 6: Choosing m' 21 and rji and then applying equations (48) and (49), we determine 
the final values of Ri and Rp. The process of computation is an iterative one and requires 
for the solution a first guess of parameters Ri and Rp. 
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6 Tooth Contact Analysis 


The purpose of TCA is to determine the influence of misalignment on the shift of the bearing 
contact and the transmission errors. This goal is to be obtained by simulation of meshing 
and contact of the pinion and gear tooth surfaces of a misaligned gear drive. 

We consider that the pinion and gear tooth surfaces are analytically represented in co- 
ordinate systems S\ and S 2 (see sections 3 and 4, respectively). The meshing of pinion and 
gear tooth surfaces is considered in fixed coordinate system Sh (figs. 11 and 12). Auxiliary 
fixed coordinate system S a and S e are applied to describe the installment of the pinion with 
respect to Sh (fig- 1 1 )• The pinion alignment error AA P is the pinion axial displacement. The 
misaligned pinion in the process of meshing with the gear performs rotation about Z c -axis. 
The current angle of rotation of the pinion is designated by <f> 1 (fig. 11). 

Auxiliary coordinate systems Sh, S c and Sj are applied to describe the installment of 
misaligned gear with respect to Sh- The errors of alignment are: the change A7 of the shaft 
angle (fig. 12), the offset A E and the gear axial displacement A A g (fig. 13). The misaligned 
gear performs rotation about the Z<f- axis, and <j> 2 is the current angle of the gear rotation. 

A TCA computer program was developed to simulate the meshing of pinion-gear tooth 
surfaces of the misaligned gear drive. The development of the TCA program is based on the 
following ideas: 

Step 1: We consider that the pinion and gear tooth surfaces and the surface unit normals 
are represented in coordinate system S\ and S 2 by vector functions 


r i(0p> V’i) and r 2 (9 g ,ip 2 ) 

(59) 

n 1 (^ P , V’i) and n 2 (0,,V>2) 

(60) 


where (9 P , t/’i) and (9 g , xp 2 ) are the surface parameters. 

Step 2: We represent now the pinion-gear tooth surfaces and their surface unit normals 
in coordinate system Sh- and take into account that the surfaces are in continuous tangency. 
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Then we obtain the following equations 


i,<M - r ! 2) (0 3 ^2>2) = 0 

-nt 2) (^s>^2,^2) = 0 


where 


rl 2) (^,^2,<A 2 ) = M H2 (<h)T 2 (9 91 ih) 

n h\ 6 P^ uti) = Lhiifa) ni(0 p ,Vh) 


n l. 2) (^9)^’2, ^ 2 ) = Lw(^2)n2(^g, ^ 2 ) 




cos <^1 sin <^1 0 0 

— sin <^1 cos <^1 0 0 

0 0 1 AA P 

0 0 0 1 


M m (&) = 


— sin A7 

0 

— cos A7 

0 ' 


' COS <£2 

— sin (f> 2 

0 

0 ' 

0 

1 

0 

A£ 


sin <f> 2 

cos <f> 2 

0 

0 

cos A7 

0 

— sin A7 

0 


0 

0 

1 

AA g 

0 

0 

0 

1 


0 

0 

0 

1 


(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 


( 68 ) 


Equations (61) and (62) represent the conditions that the contacting surfaces at the point 


of tangency have a common position vector and a common surface unit normal. Equations 
(61) and (62) yield a system of five independent scalar equations of the following structure 


/t'(0 P , = 0 /, € C 1 (i = 1..5) 


(69) 


As a reminder, vector equation (62) yields only two independent scalar equations, and 
not three, since |nj^| = |n^| = 1. 

Step 3: System (69) of five nonlinear equations contains six unknowns, but one of the 
unknowns, say <f>\, may be considered as the input parameter. Our goal is the numerical 
solution of nonlinear equations (69) by functions 


Wi), tWi), M<h)} e c 1 


(70) 
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The sought-for numerical solution is an iterative process that requires on each iteration 
the observation of the following conditions [2, 4, 5, 11]: 

(i) There is a set of parameters (the first guess) 

P(fy\li 0 \<A e \% 0 \rl>i 0) A 0) ) ( 71 ) 


that satisfies the equation system (69). 

(ii) The Jacobian taken at P differs from zero. Thus, we have 

a £ih± hi hi h) / n 


(72) 


Then, as it follows from the Theorem of Implicit Function System Existence, equation 
system (69) can be solved in the neighborhood of P by functions (70). 

Using the obtained solution, we can determine the path of contact on the pinion-gear 
tooth surface, and the transmission errors caused by misalignment. The path of contact on 
surface E, (i = 1,2) is determined by the expressions 


ri(0 P ,^i), 

e P (<f> i), 


(73) 

r 2 {0 g ,ip2), 

Wi), 

^ 2 (^ 1 ) 

(74) 

The transmission errors are determined by the equation 


A<f>2 = 

<^ 2 (^ 1 ) — 

% 

(75) 


The dimensions and orientation of the instantaneous contact ellipse at the contact point 
may be determined considering that the principal curvatures and directions of the contacting 
surfaces, and the elastic approach of the surface [2] are known (see step 7 in section 5). 
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7 Avoidance of Pinion Undercutting 


In most cases undercutting can be avoided, if the appearance of singular points on the 
generated surface is avoided. Singularities on the surface occur when the normal to the 
surface becomes equal to zero. To avoid undercutting of the pinion by the generating tool, 
the approach developed in [7, 8, 9] is applied: 

Step 1: Consider that the surface of the generating tool is represented as 

=r u (X P ,0p) (76) 


The equation of meshing is represented as 


/i(Ap,0 p ,0i) = 0 

Step 2: It is proven in [7, 8, 9] that singular points occur if 

v (‘i) q. = 0 


(77) 


(78) 


where is the velocity of the contact point in its motion over the tool surface, and vOi 1 ) 
is the relative velocity. This yields that a matrix 


A = 


dr tl dr tl 
d\ p dO v 


— v 


(*ii) 


dh df i dh drl > i 


(79) 


dX p d6 p dip i dt 

has the rank r = 2 and therefore three determinants Aj ( i = 1,2,3) of the third order must 
be equal to zero. Then we obtain that 

Fi(X p , 6 P , xpi) = Aj + Aj + = 0 (80) 

Equations (77) and (80) permit the function A P (0 P ) to be determined for the limiting 
line on the tool surface. Then, we are able to determine the limiting line on the generating 
surface by the equation 

r<,=r„(«„M«p)) < 81 > 
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Fig. 14 shows the limiting line on the pinion tool surface. 

Step 3: Using coordinate transformation, we may determine the line of singular points 
on the pinion tooth surface, (fig. 15) 

Step 4: To avoid undercutting, we have to limit the dimension of the dedendum of the 
pinion tooth. 

Fig. 16 shows the axial section of the pinion head-cutter. Parameter h represents the 
distance of a point of the axial section from the reference circle determined as 

R\ sin a p — Ri sin(ofp + A p ) = h (82) 

To verify that undercutting has been avoided the following inequality must be observed 

f?i [sin o.p - sin(a p + A P (0 P ))] > h d (83) 

where hd is the dedendum height of the pinion, and A p (0 p ) represents the function that 
corresponds to the points of the limiting line. 

The design of spiral bevel gears is based on application of special tooth element propor- 
tions for the avoidance of undercutting: small pinion dedendums and long pinion addendums. 
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8 Numerical Example 

As a numerical example the blank data is given in Table 1 . 

The gear head-cutter is a cone (figs. 2 , 3 and 4), the cutter radius is designated by R g 
(fig. 1 ), the radial setting of the head-cutter is | 0 C 2 O t2 j (fig- 5(b)), and the installment angle 
is g 2 (fig. 5). The data for the gear head-cutter that generates the gear concave side are 
presented in Table 2. 

The parameters of the pinion head-cutter were determined by application of the method 
of the local synthesis (section 5). The data for the pinion head-cutter that generates the 
pinion convex side are represented in Table 3. We considered in the numerical examples 
the meshing of the gear tooth concave side with the pinion tooth convex side. Case 1 
corresponds to the orientation of the bearing contact across the surface, case 2 corresponds 
to the orientation of the bearing contact in the longitudinal direction. 

The application of TCA for the simulation of meshing and contact permits the deter- 
mination of misalignment effects on the transmission errors and the shift of the bearing 
contact. It has been shown that in the case of application of ideal generating surfaces (with- 
out mismatch, figs. 1 and 2 ) the errors of misalignment cause indeed discontinuous almost 
linear transmission errors as shown in fig. 17 for shaft angle error A 7 . Similar functions of 
transmission errors are caused by errors A A p , A A g and A E. Table 4 shows the maximum 
transmission errors caused by misalignment. 

The results of TCA for the properly mismatched generating surfaces (see section 5) con- 
firmed that a predesigned parabolic function indeed absorbs the transmission errors caused 
by misalignment, and the resulting function of transmission is a parabolic one (fig. 18). The 
absorption of linear function of transmission errors is carried out as well in other cases of 
misalignment: A A p , A A g and A E. The bearing contact of the drive is stable, and its shift 
is permissible (fig. 19). Model 2 of the gear drive (with longitudinal direction of the bearing 
contact) is preferable due to the lower level of transmission errors caused by misalignment. 
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9 Conclusion 


From the conducted study the following general conclusions can be drawn: 

(1) An approach has been developed for the synthesis of spiral bevel gears that provides 
(i) localized bearing contact, and (ii) low level of transmission errors of a parabolic type. 
The developed approach permits two possible directions of the bearing contact: across the 
tooth surface or in the longitudinal direction. 

(2) A Tooth Contact Analysis (TCA) computer program for the investigation of the 
influence of misalignment on the shift of the bearing contact was developed. 

(3) A low level of transmission errors, the parabolic type of the function of transmission 
errors, and the localization of the bearing contact are achieved by the proper mismatch of 
contacting surfaces. 

(4) The influence of the following errors of alignment was investigated: (i) for axial 
displacement of the pinion, (ii) axial displacement of the gear, (iii) offset, and (iv) change 
of the shaft angle. These types of misalignment were proven to cause discontinuous almost 
linear functions of transmission errors, but they are absorbed by the predesigned parabolic 
function of transmission errors. 

(5) Conditions of nonundercutting of the pinion were determined. 

The results of this investigation show that a predesigned parabolic function can indeed 
absorb the linear functions of transmission errors caused by misalignment. The design of 
gears with a longitudinal bearing contact (in comparison with the bearing contact across the 
surface) is preferable since a lower level of transmission errors can be obtained. 
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10 Manual for Computer Program 


Program Names and Purpose of the Programs 
There are three programs 

1. Program for Local Synthesis: Localsyn.for 

2. Program for TCA: Tca.for 

3. Program for Undercutting: Undercut.for 

These programs are directed at the synthesis of the spiral bevel gear with uniform tooth 
height by using mismatched generating surfaces. The programs cover the local synthesis, 
tooth contact analysis and nonundercutting conditions. Using the programs, one can obtain 
the tooth surfaces, the contact lines on the tooth surface, the contact path on the tooth 
surface, the transmission errors and the bearing contact caused by misalignment of the gear 
drive, and the limiting lines on the generating tool surface and the pinion tooth surface. 

Environment for Running the Programs 

These programs were developed by application on an IBM PC and can be run using the 
software ” Power Fortran”. 

An application of the subroutine HYBRD1 [11] for solving a system of nonlinear equations 
and several other subroutines that was called by HYBRD1 are required and included. 

Input Data 
1. Blank data 


TNI — Pinion number of teeth 
TN2 — Gear number of teeth 
TW — Face width of gear (mm) 
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GAMA — Shaft angle (degree) 

Betal — Pinion spiral angle (degree) 

Beta2 — Gear spiral angle (degree) 

EllipseDelta — Elastic approach (mm) 

2. Gear cutter specification 

RU2 — Gear nominal cutter radius (mm) 

PW2 — Point width of gear cutter (mm) 

AFA-g — Blade angle of gear cutter (degree) 

Rg — Cutter radius (mm) 

3. Gear machine-tool settings 

GAMA2 — Gear machine pitch angle (degree) 
Sr2 — Radial setting (mm) 
q2 — Installment angle (degree) 

4. Pinion machine-tool settings 

GAMA1 — Pinion machine pitch angle (degree) 
Sri — Radial setting (mm) 
ql — Installment angle (degree) 

5. Pinion cutter specification 
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Rp — Cutter radius (mm) 

R1 — Radius of surface of revolution(mm) 

AFA 4 > — Profile angle of gear cutter (degree) 

6. Misalignments 

H — Axial displacement of the pinion (mm) 

Q — Axial displacement of the gear (mm) 

V — Offset displacement (mm) 

Delta — Change of shaft angle (arc min.) 

7. Local synthesis 

Etal — Tangent to the contact path on pinion surface at the mean contact point(degree) 
plantm21 — Coefficient of the parabolic function 

Output data files 

File philphi2.kl: Transmission errors Afo 
File sgthetag.kl: Contact line on gear generating surface 
File spthetap.kl: Contact line on pinion generating surface 
File contactp.kl: Contact path on pinion tooth surface 
File contactg.kl: Contact path on gear tooth surface 
File ellipse. kl: Bearing contact on the pinion surface 

File undercut. kl: The limiting line on pinion tooth surface (for avoidanceof undercutting) 
Procedure of using the programs 

Step 1: Run program Tea. for for the condition of no misalignment by supplying the first 
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guess of R1 and Rp. 

Step 2: In the output file contactg.&f we will get from the first line the x, y and z 
coordinates of the first contact point. 

Step 3: Run program Localsyn, for for the desired plantm21 and etal. Then we can get 
R1 and Rp. 

Step 4: Check if R1 and Rp at Step 1 and Step 4 are the same or not. If both are the 
same then go to Step 6. 

Step 5: Use the new values that we got from Step 4, recalculate the first contact point 
by running program Tea. for, and go to Step 3. 

Step 6: Run program Tea. for with misalignment to obtain the transmission errors in the 
output file philphi2.&/. 

Step 7: Run the program Undercut. for to check up the undercutting in the output file 
underp.fc/. 
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Table 1 : Blank Data 



Pinion 

Gear 

Ni,N 2 , Number of teeth 

11 

41 

7 , Shaft angle 

O 

O 

CT> 


Mean spiral angle 

o 

lO 

CO 

35° 

Hand of spiral 


LH 

Whole depth (mm) 

6.5 

6.5 

Tooth module (mm) 

4.33 


Face width (mm) 

27.25 

27.25 

7 i, 72 , Pitch angles 




Table 2: Parameters and Installment of Gear Head-Cutter on gear concave side 


a g , Blade angle 

O 

O 

R g , Cutter radius at mean point (mm) 

78.52 

S t2 , Radial setting (mm) 

70.53 

q 2 , Installment angle 

-62° 14' 


Table 3: Parameters and Installment of the Pinion Head-Cutter on pinion convex side 



Case 1 

Case 2 

a p , Profile angle 

20 ° 

20 ° 

INPUT 



rji, Tangent direction of the contact path 

171° 

92° 

Derivative of 4>2{<t>\) 

-1.3e-3 

-1.2e-3 

A^ 2 = 0 . 5 m 2 1 ( 7 r/A f i) 2 , Theoretical Max. (”) 

-10.94 

-10.09 

OUTPUT 



M, Mean contact point in S m (mm) 

(79.88, 0.39, 0.17) 

(77.83, 1.64, 0.72) 

Up, Cutter radius at mean point (mm) 

78.0 

64.7 

Ri, Radius of the surface of revolution (mm) 

235.0 

765.0 

Sri (mm) 

70.30 

65.38 

< 7 !, Installment angle 

-6i o 5r 

— 51°24' 

Length of major axis of contact ellipse (mm) 

12.54 

4.5 
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Table 4: Maximum Transmission Errors for Generating Surfaces with Mismatch 



Afo in arc sec . 


Case 1 



8.8 

16.2 


11.5 

12.5 

A E = 0.1mm 

11 

15 

A 7 = 3' 

10.7 

13.5 
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Contact Path 




Fig. 2: Generating cone and generating surface of revolution 
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Fig. 3: Mismatched generating surfaces 


Generating 

cones 
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Gear pitch cone 






Fig. 7: (a)Convex (Inside blade) and (b)concave (outside blade) sides of the generating 
blades and generating surfaces of revolution 
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Fig. 8: Coordinate systems S m , S a and Si 


(a) 


(b) 



Fig. 9: Transmission function and predesigned parabolic function of transmission errors, 
-pinion rotation angle; <j> 2 -gear rotation angle; A ^-transmission error 
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Fig. 15: Limiting line on pinion tooth surface 



(Arc sec) 


z *.1 



Fig. 16: For the derivation of the limiting value of the dedendum 



(j) 1 (Deg) 

Fig. 17: Transmission errors for a misaligned gear drive with ideal surfaces: A7 = 3 arc 
min. 
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